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ABSTRACT 


Future solar sail missions will require sails with dimen- 
sions on the order of 100 m to 1 km. At these sizes, given 
the gossamer nature of the sail supporting structures, flex- 
ible modes may be low enough to interact with the control 
system. This paper develops a practical analysis of the 
flexible interactions using state-space systems and modal 
data from standard finite element models of the sail sub- 
system. The modal data is combined with a rigid core bus 
to create a modal coordinate state-space plant, which can 
be analyzed for stability with a state-space controller. Re- 
sults are presented for an 80 m sail for both collocated ac- 
tuation and control by actuators mounted at the sail tips. 


1. INTRODUCTION 


Future solar sail missions, such as NASA’s Solar Polar 
Imager, will require sails with dimensions on the order of 
100 m to 1 km. At these sizes, given the gossamer na- 
ture of the sail supporting structures, flexible modes may 
be low enough to interact with the control system. This 
paper develops a practical analysis of the flexible inter- 
actions using state-space systems and modal data from 
finite element models of the sail. The tools, part of the 
MATLAB™-based Solar Sail Module for the Spacecraft 
Control Toolbox developed under the NASA In-Space 
Propulsion program, include a graphical user interface 
which allows users to quickly perform analyses (Thomas 
et al. 2004). 

Standard finite-element models of the flexible sail body 
are loaded and the modal data is used to create a modal 
coordinate state-space system. The analyst specifies 
which modes to include, which nodes are of interest for 
force inputs and displacement outputs, any nominal mo- 
mentum in the system, and any steady rates. The system 
is linearized about the nominal attitude and body rates. 
The state-space plant can then be analyzed with a state- 
space controller, and Bode, Nyquist, step and impulse re- 
sponses generated. The approach is general for any rigid 


core with a flexible appendage. 

The sail analyzed is a square sail with four gossamer 
masts and four triangular sail membranes suspended be- 
tween them. The coordinate frame of the sail - shown in 
a sun-synchronous Earth orbit - is given in Figure 1 . 



Figure 1. Sail Coordinate System (ATK Space Systems) 


A finite element model of the scalable sail has been pro- 
vided for this analysis by ATK Space Systems (Goleta, 
CA) (Murphy et al. 2005). The first 30 modes range from 
about 0.04 to 0.08 Hz. The model and modes can be ef- 
fectively visualized in MATLAB as seen in Figure 2 on 
the next page. One control method under consideration 
is a suite of pulsed plasma thrusters mounted at the tip 
of the masts (Wie et al. 2004). On the +F axis for this 
model, the corresponding node at which the force will be 
applied is 50041. A +X force applied here will produce 
a negative rotation about the Z axis. 

This paper presents results for an 80 m sail controlled by 
a PID, such as would be implemented by reaction wheels 
or control moment gyros (CMGs) on the core bus. In ad- 
dition the response to actuators mounted on the sail tips is 
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Figure 2 . Solar sail mode visualization in MATLAB 


shown, such as by control vanes or small thrusters. Both 
the overall control stability and the resulting deflection of 
any node is obtained. 


2. DERIVATION OF LINEAR PLANT MODEL 


We require a state-space form of the dynamics, i.e. 

x — Ax 4- Bu 
y = Cx 

where x is the state, u is the control input, and y is the 
measured output. 

Our vehicle consists of a rigid core and a flexible ap- 
pendage. We begin by defining a set of nodes on the 
appendage. If the vector from system center of mass 
(C sys ) to the origin is c, the nominal node locations are 
rk, and the nodal displacements are p*, then the vectors 
from C sys to the nodes are dk = c 4- r*. + p fc . This is 
depicted in Figure 3. 

We also define the modal displacements 77, where the 
transformation matrix $ gives the nodal displacement per 
unit modal displacement for each mode. $ is constructed 
as a 3 n x m matrix, where n is the number of nodes and 
m is the number of modes. This can be written compactly 
in matrix notation as 

P = $T} ( 1 ) 

Our goal is to develop state-space equations using the 
modal coordinates 77. We begin with the total angular 
momentum of a rigid core with a flexible body about the 
spacecraft center-of-mass. Note that the sum includes the 
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Figure 3 . Node vector definitions 


displacement of the core from the system center of mass 
as the sum starts at 0. 

n 

H = A(Icj) + Y j ™ k D£ D k ( 2 ) 

*=0 

where A is the rotation from the body frame to the inertial 
frame, I is the inertia matrix of the core and Dk = Adk 
are the node vectors from the system center of mass in the 
inertial frame. Taking the inertial derivative and trans- 
forming to the body frame we have 

n 

T = Id) -f u ) x Iu) 4 * y irikd £ dk ( 3 ) 

k-0 

where T is the external torque on the spacecraft. 

The second key equation is the inertial equation of motion 
for the nodal masses, including the motion of the space- 
craft center-of-mass, X , produced by the total external 
force. 

fk = mk(dk+X) ( 4 ) 

where fk is the sum of the internal forces exerted by ad- 
jacent nodes and any external forces on the nodes and 

We manipulate this second equation first. The definition 
and derivatives of d are 

d = r -f c-f- p 
d = u; x d-hp'fc 

d = t j x u x d F d) x d 2c 0 x (p + c) -f p -f c 

where c, taking the origin O at the core center so that ro 
is zero, is 

c = (5) 

' rriT 

/c=l 1 
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Substituting this into Equation 4, dropping the nonlinear 
terms in u, p and c, 

Y F 

fk W u + Pk + C + = — ) (6) 

2 v m 

The internal mass forces can also be represented as a 
spring force, giving 

/* = ~~ y: kkj (pj — Pfc) + Fk = —Kp 4- F* (7) 
i 

where A: is a physical property of the structure and Fk is 
the external force on each node. 

Substituting the definitions of fk and c we have 


where Z is the modal damping matrix and u are the nodal 
inputs, except that having M will result in the combined 
model free frequencies being higher than the appendage 
cantilevered frequencies. 

Next we move to development of the torque equation (3). 
First we write an intermediate step of the equation by sub- 
stituting in the equation for d, neglecting only the nonlin- 
ear terms in p and c: 

T = Ilj + cj*Iuj (13) 

n 

- (u^u^dk -d%uj + Pk 4-c) 

k~0 

The total inertia It is defined as 


mk(pk+c) + Kp = 


V F 

~rrik(—d^u) -f ~ — ) 4- Fk (8) 

2^771 


-rrik 


+ F k = 


rriT - rrifc 


^ Tn rriT 

= v k F k - n k Yl Fj 


F k -m k \^ — 


pk 4- c 


fo-Y 


Pj 

rriT 


771 J 1 — 771 

m>T 


h-Y, 


m 3 Pj 
rriT 


- vkh-'FjVjPj 

m k {v k p k - J2PjPi) + Kp = 
m k d£tjj + u k F k - Mfc F i 


(9) 


where K is symmetric and full. For three-dimensional 
displacements, p is 3nx 1 and K is 3nx 3n. Now we 
assign M to be the coefficients of p and N to be the co- 
efficients of F, and we have in matrix form 

Mp 4- Kp = Md*w 4- NF (10) 


$ for the cantilevered flex modes is normalized so that 
the modal masses are equal to one and 

= E 


= n 2 

Therefore, switching now to modal coordinates and mak- 
ing the above substitutions in Equation 10 we have 


M&j + K9ri = Md x v + NF 
<f> T M$Ti + $ T K$T] = $ r (M d* CJ + NF) 
$ T M$ii + n 2 v = 9 T (Md*u + NF)(ll) 

This is comparable to the general form of a modal dy- 
namics equation for a flexible system (including damp- 
ing) (Gawronski 1998), 

77 4- 2Zfb) 4 - Q 2 r? = Bu (12) 


n 

It ~ I core ^ TTLkd^ 

k-0 


We gather terms which contribute to Jr, note that the sum 
of the terms involving c is zero and p 0 is also zero, and 
convert to modal coordinates: 

n 

T = Itv + (j*Itw - '^mkd£$fj (14) 

fc=i 

The torque T is the sum of all external torques, including 
the torque onthe core T c and any forces applied at the 
nodes. 

T = T c + Y^dZF k (15) 

We can now write both the core and mass equations in 
matrix form with the terms containing u and 77 on the 
left and the included modes 1-m listed explicitly. Let 
Fm = Zk m kdk <t?k* ^ere <fP k is the 3 x 1 displacement 
vector for node k due to mode j. These terms indicate 
the contribution for each mode, which is a sum over the 
contributions from each node. 
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Grouping the modal contributions into one matrix I m and 
explicitly writing out the core and node inputs we have 
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(17) 


To obtain a state-space form from this equation, we will 
have to divide by the general inertia matrix on the left, 
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I g . The output matrix C is composed of the submatrices 
$ for the output nodes of interest. 


A 

B 
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I S \ 


I 9 \ 
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0 

-n 2 


e 3 
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0 0 
0 m 0 


0 

-2ZO 

« ‘ 
$ T 


(18) 

(19) 

( 20 ) 


In addition, we have left in the (Icu) x nonlinear term 
for now, which will need to be linearized. If all rates 
are small, the term can be neglected. Otherwise, there 
are two common cases: first the presence of a large bias 
momentum (such as a momentum wheel), and second a 
steady body rate. In the case of bias momentum and small 
body rates, the bias drives the nonlinear term and the lin- 
earized version is simply 


CJ X Iu « —hgLU 

In the case of a steady rate ujq, we linearize around the 
rate and drop terms of second order. 


v x Iv = (a; -f cjq) x I(u + cjo) 

~ (—(I u>o) x 4* Vq I)gj -f Vq /u>o 


To complete our state-space representation we need to 
add states for the attitude kinematics. We consider three 
easily linearized cases, small Euler angles, small quater- 
nion deviations, and perturbations from a rotating frame. 
Small 3-2-1 Euler angles can be integrated to first order 
directly from the body rates, giving 

e = [ o E] \ 9 

L J V 

If the spacecraft is undergoing a steady rotation and we 
designate the nominal rate vector as cu 0 , then the angular 
velocity is simplified with the small angle approximation 
as follows, 

v « 0 + (E - 0 x )u Q 
Rearranging to solve for 0 we have 
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For small angle quaternions, the first component ( q s ) is 
constant and the derivative reduces to one-half the body 
rates for the latter three terms (q r \ 
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3. SIMPLE NON-COLLOCATED SYSTEMS 


When a control vane or thruster is used at the tip of the 
solar sail for control the actuator is separated from the 
sensor by a flexible structure. If the bandwidth of the con- 
troller is much lower than the lowest flexible frequencies, 
as is the case with the solar sail control system discussed 
in this paper, this does not pose a problem. If the band- 
width approaches that of the bending modes then we have 
the classic noi^follocated actuator/sensor problem. This 
section provides an introduction to the theory of control 
of systems with noncollocated actuators and sensors, an 
area that has been studied extensively in the context of 
large space structures control. 


3.1. One-dimensional system 


We begin our control analysis by compensating a simple 
1 dimensional two-mass, one-spring system with a non- 
collocated sensor and actuator. The spring has a stiffness 
k and the system has damping c. 


*1. Fi x 2 , F 2 
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m x x i + c(x i - x 2 ) 4- k(x i - x 2 ) = F l (21) 
m 2 x 2 + c(x 2 -ii) + k{x 2 -xi) = F 2 


In the frequency domain, neglecting damping for the mo- 
ment, the system becomes 




s 2 + w 2 w\ 

Fi/mi 

rci ' 


U)f s 2 + w\ 

Film-z 

X2 J 

s 2 (s 2 -1- (u; 2 + iu|)) 


where v\ = k/mi and = k/m 2 . From this we see 
that for the collocated case F 2 to x 2 , at low and high fre- 
quencies relative to the flex resonance the system acts like 
a rigid body, 1/s 2 . The lead from the zero(s) is cancelled 
by the lag from the pole(s) so that the phase after the flex 
frequencies remains 180°. For F 2 to x lf however, at high 
frequencies the denominator approximates 1/s 4 , result- 
ing in permanent phase lag. This is easily seen in the 
open loop response for both cases in Figure 4 on the next 
page. 

The system properties are selected to give a natural fre- 
quency near the lowest sail mode, mi is 10 kg, m 2 is 
1 kg, k is 0.3, c is 0.03, v n is 0.574 rad/s and $> = 
[0.0953, -0.953]. 
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Figure 4. Open loop response for collocated (left) and non-collocated (right) actuation , ID system 


The state is [x v] T . The state space matrices are simple. 
First we define the stiffness, damping and mass matrices 
K, D and M, 


-0.021124 

-0.0061566 

0 



There is now a tradeoff between amount of phase lead and 
excitation of the flex mode. Figure 5 on the following 
page shows the open loop and step responses for phase 
leads of 40° and 65°. 


M 


771 1 0 

0 7712 


and then we can write out the state space. 


A 

B 

C 


0 E 

- M —M~ l D 

' 0 
M - 1 

[E 0 ] 


We first apply a forward gain for the desired transient re- 
sponse, K = 0.01. The result is unstable as shown by the 
eigenvalues, 


-0.00165 + 
-0.00165 - 
1 . 2661e-08 + 
1 . 2661e-08 - 
0 


0 . 57366i 
0 . 57366i 
0 . 030193i 
0 . 030193i 


Next we add a notch filter to reduce the gain seen by the 
flex mode. This should allow us to use the higher phase 
lead of 65° (i.e. more damping) with less excitation of the 
mode. We set the notch frequency equal to the resonant 
frequency (0.574 rad/s) with a gain reduction of 20 dB 
and a notch width of 0.1 rad/s. See Figure 6 on the next 
page, and note that the step response is lacking the small 
oscillation seen in Figure 5 on the following page. 

Lastly, we need to add an integrator to drive the offset to 
zero. We will need this in the solar sail controller since 
the pointing during delta-V changes controls the thrust 
vector; we essentially want to limit the average pointing 
error over a maneuver within some tolerance. There is 
a trade between the time constant of the integrator and 
the forward gain; an infinitely high gain would result in 
a negligible offset, and correspondingly adding the inte- 
grator should allow for a higher gain. We will use a slow 
integrator of 3600 s, and add it in parallel with the lead 
compensator and filter. See Figure 7 on page 7 for the 
step response before and after the addition. 


We then add phase lead at the crossover frequency for the 
desired damping, with a maximum phase shift of 40° at 
Wmax = 1.0~ 2 . This pushes the crossover frequency up, 
so we lower the gain to 0.001 and raise w m ax to 2.0 , 
which gives a stable system with the eigenvalues 


This same loop shaping procedure can be used with more 
complex systems such as the solar sail. 


-0.0014114 + 0 . 57287i 

-0.0014114 - 0 . 57287i 

-0.063411 




Magnitude (db) 












Figure 7. Step response of ID system before and after integral is added . 


3.2. Planar system 


Now that we have established the response and control of 
the simplest flex system, we check our theoretical devel- 
opment for n flex nodes with a rotational planar system 
consisting of a core and two masses, having two flex fre- 
quencies. The core mass is 100 kg, mi is 2 kg and ri is 
2 m, and m 2 is 1 kg and r 2 is 4 m, and the flex displace- 
ments are transverse to the node locations, as shown in 
Figure 8. For a spring stiffness of 0.5, the cantilevered 
natural frequencies are 0.38268 and 0.92388 rad/s. If the 
flex body is correctly linked to the core, we expect that 
the response to a force at r 2 will be 4 times the response 
to a unit torque on the core. We see in Figure 9 on the 
next page that we do get the collocated and noncollocated 
responses we expect. Note the upward shift of the flex 
frequencies from the cantilevered frequencies to the free 
frequencies. 
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Figure 8. Simple rotational system 


4. SOLAR SAIL ANALYSIS 


In the case of the 80 m sail, the flex frequencies are rel- 
atively high compared to our desired control bandwidth 


u> n of 0.002 rad/s(Wie et al. 2005), but also very close 
together, with the first 30 modes between 0.25 and 0.5 
rad/s. The many overlapping poles and zeros may pre- 
vent us from seeing the clear noncollocated response as 
in the simple models. 

First we check the open-loop sail plant as developed in 
Section 2, starting with the collocated response. Next 


Input: Tz to Output: thz 




Figure 10, Sail plant response with modes 1-30 


we consider our desired force input at a mast tip. We 
choose a node (50041) on the +Y boom and consider 
forces in the +X and -f Z directions, which will produce 
respectively ~T Z and +T S . The inertia about the X axis 
is twice that of the other axes, and the boom tip is about 
56 m from the core center. This is a complex system, 
and we do not see a simple noncollocated response in 
Figure 1 1 on the following page. We therefore also look 
at a node in the middle of a sail membrane (598). When 
including modes 1-4, we get a phase drop of 180° and a 
rigid body torque response (Figure 12 on page 9). Note 
the nodal output in the step response plots. 
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figure 9. Response for a core with a 2 node flex body (planar 
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Figure 11. Sail plant response of sail to mast tip 
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actuation (node 50041) and membrane actuation (node 598). including 
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Figure 12. Membrane actuation ( node 598) with PID controller, u) n = 0.001 rad/s 


We are interested in analyzing controllers at different 
bandwidths using our modal plant model. We first ap- 
ply a PID controller, with modes 1-30 in the plant, and 
check the eigenvalues for stability. Both collocated ac- 
tuation (core torques to core angles) and actuation at the 
mast tips are stable for all u; n tried from 0.001 to 1 rad/s 
(with a damping ratio C of 0.8 and integrator time con- 
stant of 3600 s. However when the node in the membrane 
is considered, the system becomes unstable above 0.002 
rad/s, so this provides an opportunity to try the loop shap- 
ing procedure outlined in Section 3 on page 4. Figure 12 
shows the response with a PID controller at 0.001 rad/s 
for comparison. Note the crossover at about 0.02 rad/s. 

Following the same procedure, we use a gain of 0.6, and 
a lead compensator with a maximum phase frequency of 
0.12 rad/s and phase gain of 60°. Then we add a notch 
filter at the apparent frequency where the modes are clus- 
tered, about 1.5 rad/s, with a half-notch width of 1 rad/s 
and a gain drop of -30 dB. The frequency responses dur- 
ing this process and resulting stable step response are 
shown in Figure 13 on the following page. Note that we 
this compensator has a crossover frequency of about 0.1 
rad/s. 

Analyzing this sail configuration reminds us of an im- 
portant fact, that the noncollocated response is irrele- 
vant if the system is very stiff. In this case, the masts 
are very stiff relative to the sail membranes. Within the 
first 30 modes, only several have any boom involvement 
and those have very small displacement magnitudes. Al- 
though we might expect stability problems when per- 
forming actuation at the end of very long masts - here, 
56 m - in fact there are none. However, actuating nodes 
on the membranes provides the expected response, both 
validating the analysis and indicating that if attitude sens- 
ing were for some reason placed on the membranes, or if 
in other configurations the supporting structure is equally 
flexible relative to the membrane modes, case must be 
taken in designing the control law. 
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